Radial Velocity Jitter in Stars from the California and Carnegie 
Planet Search at Keck Observatory 1 
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ABSTRACT 

I present an empirical model for predicting a star's radial velocity jitter from its B — V color, 
activity level, and absolute magnitude. This model is based on observations of ~ 450 well- 
observed stars from Keck Observatory for the California and Carnegie Planet Search Program. 
The model includes noise from both astrophysical sources and systematic errors, and describes 
jitter as generally increasing with a star's activity and height above the main sequence. 

Subject headings: stars: activity, techniques: radial velocities 



1. Introduction 

The California and Carnegie Planet Search has 
been observing ~ 1200 stars, some since 1997, at 
Keck Observatory in an ongoing program to de- 
tect radial velocity variations caused by planetary 
companions. One source of error in the measured 
velocities is "jitter" due, in part, to flows and in- 
homogeneities on the stellar surface that can pro- 
duce variations in the measured radial velocity of 
a star, and can even mimic the signature of plan- 
etary companions (Henry, Donahue, & Baliunas 
(2002); Queloz et al. (2001); Santos et al. (2003)). 
Marcy et al. (2005) describe how the California 
and Carnegie Planet Search distinguishes plane- 
tary signals from jitter. 

The best-fit orbital parameters of a Keplerian 
velocity model depend on the relative errors of 
points in an observed radial velocity time series. 
Thus, one requires a quantitative model of all 
sources of noise, including jitter, in order to prop- 
erly calculate the orbits of a potential planetary 
companion. 

Saar & Donahue (1997) studied jitter in F, G, 



1 Based on observations obtained at the W. M. Keck 
Observatory, which is operated jointly by the University 
of California and the California Institute of Technology. 
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and K dwarfs, modeling the effects of spots and 
convection on radial velocity measurements as a 
function of T c g and activity. They concluded that 
radial velocity jitter should increase with T g, ac- 
tivity, and v sin i, and found this to be the case em- 
pirically in a small sample of active stars. Paulson, 
Cochran, & Hatzes (2004) found that among four 
Hyades stars with significant radial velocity vari- 
ation, three showed correlated photometric vari- 
ation, consistent with the effects of rotationally 
modulated spots and plage. Saar, Butler, & Marcy 
(1998) used a sample from the Lick Planet Search 
results to produce a formula for predicting jitter 
in stars given -P ro t, usini, i>mac, and photometric 
variability. 

The purpose of this work is not to explore the 
astrophysical sources of jitter or to refine the for- 
mula of Saar, Butler, & Marcy (1998) which pre- 
dicts jitter from photometric variability and v sini 
(a quantity that is difficult to measure in stars as 
old as typical Planet Search stars), but to pro- 
vide an empirical estimate of jitter using stellar 
parameters available without an independent pho- 
tometric program or detailed spectroscopic analy- 
sis: My, B — V, and S. The strategy is to divide 
the Planet Search sample into groups of stars with 
similar values for these parameters. The model 
then simply predicts that a given star will exhibit 
levels of jitter similar to those stars with similar 
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properties, i.e., in the same group. 

2. Data 

2.1. The Sample 

To measure the radial velocity jitter this work 
employs a sample of stars that have useful radial- 
velocity time series and no detected companions. 
From the ~ 1200 stars in the Keck Planet Search 
program, this sample comprises 531 stars with at 
least nine radial velocity measurements, with ac- 
tivity levels reported in Wright et al. (2004), and 
Hipparcos colors and parallaxes (ESA 1997). 

Rather than attempt to subtract the velocity 
signatures of known companions from the sample, 
some of which may be rather poorly determined, 
the sample here excludes those stars harboring 
published or announced companions from the sam- 
ple altogether. Furthermore, this sample excludes 
stars that show clear, unambiguous evidence for 
unannounced companions (stellar and planetary). 
These cases involve very long period companions 
of unknown mass where the radial velocities show 
curvature, or Kcplerian amplitudes of more than 
1 kms -1 . Finally, the sample excludes several 
"borderline" stars whose radial-velocity time se- 
ries have best-fit Keplerians with false alarm prob- 
abilities (Marcy et al. 2005) of less than 0.1. 

The final sample comprises 448 stars. 

2.2. The Evolution Metric 

Following Wright (2004), this work employs the 
evolution metric 

AM Y = M Vt Ms(B-V)-M v (1) 

as a quantitative proxy for evolution. Here, 
My- and B — V- values are from the Hipparcos cat- 
alog, and 

Mv M 8{B-V)=Y,a i (B-V) i (2) 

represents the Hipparcos average main se- 
quence where a = {1.11255, 5.79062, -16.76829, 
76.47777,-140.08488, 127.38044, -49.71805, -8.24265, 
14.07945, -3.43155}. Fig. 1 shows M v ,ms plotted 
on a color-magnitude diagram with the stars from 
the sample described in 2.1. Evolved stars, lying 
above the main sequence, have positive values of 
AMy, while subdwarfs have negative values. 



While a more appropriate metric might employ 
measured abundances to lift the degeneracy be- 
tween the effects of metallicity and evolution on 
AMy, the goal here is to provide a description 
of radial velocity jitter without the benefit such 
measurements. 

2.3. The Activity Metric 

Active stars, especially active F dwarfs, exhibit 
higher levels of jitter than their inactive counter- 
parts. To quantify this relationship, this work em- 
ploys the activity metric S, as measured in Wright 
et al. (2004) for the Planet Search stars, which rep- 
resents the amount of emission in the Ca II H and 
K line core normalized by nearby continuum re- 
gions. Wright et al. (2004), following Noyes et al. 
(1984), transformed S to -Rhk which represents 
the fraction of a star's total luminosity emitted by 
the chromosphere in the Ca II H and K lines. Ex- 
ample spectra containing Ca II H and K lines can 
be found in most of the planet discovery papers 
by the California and Carnegie Planet Search, e.g. 
Butler et al. (2003). 

The quantity R' KK , however, is not well cali- 
brated for stars with B-V < 0.4 or B-V > 0.9. 
To overcome this, this work employs an alternative 
metric of stellar activity, F Ca jj, defined by Rutten 
(1984) (who used the expression "Fn + fk" in that 
work) and described again in Rutten (1986). This 
metric differs from R' KK primarily in that it does 
not correct for a photospheric component to the 
Ca II H and K lines and in that it is a flux density, 
rather than a fraction, since it is constructed using 
an additional factor of T c 4 ff . Rutten & Schrijver 
(1987) argued that F c jj is a more appropriate 
metric for comparing activity levels among stars 
since it is related to soft X-ray and Mg II h and k 
fluxes by power laws with no color or luminosity 
dependence. Hall & Lockwood (1995) determined 
the units of F CaII to be 7.97 x 10 5 erg cm~ 2 s _1 . 
Rutten (1984) prescribes 

^CaII = ^cfT o 4 ff 10- 14 (3) 

to transform the ^-values of the stars in the 
sample here to F Ca tt where 

logC c f = 0.25(B-Vf-l.33(B-V) 2 +0A3(B-V)+0.24: 

(4) 
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for main sequence stars with 0.3 < B — V < 1.6. 
This work takes T e s for sample stars from Valenti 
& Fischer (2005) or, for the few stars here without 
such measurements, from the T e s-B — ^relation of 
Flower (1996). 

Fig. 2 shows F Ca jj as a function of color for 
the sample stars here. Rutten (1987) tabulated 
the empirical minimum activity level F Ca tt min ob- 
served for stars as a function of B — V. Table 1 and 
the solid line in Fig. 2 give this tabulation, con- 
verted to the units of this work. This minimum 
level ostensibly represents radiative flux not re- 
lated to magnetic activity. A comparison of Fig. 2 
with Fig. 6a of Rutten (1984) shows good agree- 
ment (note that the curves in Rutten's figure have 
been superseded by the curves referenced above). 

This work employs the activity metric AF Ca jj, 
constructed as: 

AF CaII = F CaII - F CaII,min ( 5 ) 

For some metal-poor, blue stars, AF Ca tt is neg- 
ative. Apparently, these stars' very low metal- 
licities ([Fe/H] < —1) require either a different 
calibration of F Ca jj or a different expression for 
F Ca u min (or both) . This is not a problem in this 
analysis, because, as described below, AF Ca jj is 
not needed to predict the jitter for stars as inactive 
as these. 

Henceforth, this work uses the term "active" 
to refer to stars with F Ca jj > 0.6. This limit 
corresponds to logi?^ K ~ —4.95 at B — V = 0.5 
and logi?^ K 4.8 at B — V = 0.9. 

2.4. Radial Velocity Measurements 

2.4.1. The Jitter Metric 

Many sample stars exhibit linear trends in their 
radial velocity time series due to long-period plan- 
etary or stellar companions. To remove such ef- 
fects, the standard deviation of the velocities, a rv 
is defined about a linear trend, v r = a + bt, for all 
of the time series: 

2 K,i - a- bU) 2 
CT ™ = .L j^— 2 (6) 

i=0 

The quantity a rv includes random measure- 
ment uncertainties such as those arising from pho- 
ton statistics, estimated as follows. Each observa- 



Table 1: The lower boundary of the surface flux 
F Ca 11 as a function of B — V, adapted from Rut- 
ten (1987) to the units of this work. 



B — V 


Ca 11, min 


0.4 


5.24 


0.45 


3.80 


0.5 


2.88 


0.55 


2.29 


0.6 


1.82 


0.65 


1.41 


0.7 


1.09 


0.75 


0.831 


0.8 


0.645 


0.85 


0.489 


0.9 


0.363 


0.95 


0.281 


1.0 


0.218 


1.05 


0.174 


1.1 


0.135 


1.15 


0.105 


1.2 


0.079 


1.25 


0.063 


1.3 


0.052 


1.35 


0.044 


1.4 


0.038 


1.5 


0.029 


1.6 


0.025 
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tion of a star consists of a spectrum divided into 
400 chunks, each of which is analyzed indepen- 
dently the Planet Search Doppler analysis pipeline 
to determine the value of Av, the differential ra- 
dial velocity of the star with respect to a template 
observation. The distribution of the Av values of 
these 400 chunks provides an estimate of the mea- 
surement error for a given observation. Because 
the Planet Search employs an exposure meter at 
Keck, most observations of a given star have the 
same signal-to-noise (S/N) ratio, so the measure- 
ment errors for all observations of a given star are 
usually very similar. Thus, a single number er m0 as 
provides an estimate of the contribution of mea- 
surement error to u rv for observations of a given 
star. 

Following Saar, Butler, & Marcy (1998), this 
work employs the metric a' rv — (o- 2 rv — crj^J 1 / 2 
to quantify radial velocity variations in excess 
of cr moas . The exposure meter at Keck ensures 
that those velocity errors limited by S/N are 
roughly constant across the sample, contributing 
2-4 ms _1 to the error budget. The exceptions to 
this rule are faint M dwarfs, which receive a max- 
imum of ~ 10 minutes per exposure and therefore 
have photon-limited errors of 3-8 ms _1 , and ac- 
tive F dwarfs that have high astrophysical jitter 
and receive short exposure times to keep errors 
photon-limited. For these stars, the observed rms 
velocities are nearly equal to <r mcas , indicating that 
photon statistics dominate the error budget. In 
cases in which where radial velocity variations in 
excess of a meas remain undetected, this work sim- 
ply assigns an upper limit of cr moas /2: 



iV rv mcas u rv ^ 4 mcas 

5 2 
< Cmcas/2 o rv < jCT meas 

2.4-2. The Components of o~' rv 

In addition to astrophysical jitter, o~' rv also mea- 
sures undetected planetary companions and sys- 
tematic errors as discussed below. 

Many stars in this sample undoubtedly harbor 
still-undetected, low-amplitude planets. Only af- 
ter many planetary orbits with good phase cov- 
erage can one use the periodic nature of the sig- 
nals and Keplerian fits to distinguish these signals 
from jitter (Marcy et al. 2005). Refining these 



fits to determine the precise orbital elements of a 
planetary companion requires still more observa- 
tions, and the statistics of the residuals to these 
fits is complex. It is for the sake of simplicity and 
uniformity in the sample that this work employs 
only linear fits to the velocity time series and not 
Keplerians. 

Systematic errors have many sources, including 
imperfect deconvolution of template spectra (es- 
pecially for M dwarfs) and poor charge transfer in 
the HIRES Tck CCD. The Planet Search estimates 
these errors to be on the order of 3 ms _1 for M 
dwarfs and < 2 ms~ 1 for other stars (G. W. Marcy 
2005, private communication). Some component 
of these systematic errors will be accounted for 
which the subtraction of cr moas from o~ rv , but sys- 
tematic errors, by their nature, are not well char- 
acterized by Gaussian noise, so some component 
of these errors probably remain in a' rv . 

3. Analysis 

The purpose of this paper is to quantify the 
known dependencies of radial velocity jitter on 
stellar parameters, namely that active, blue, and 
evolved stars exhibit more jitter than old, main 
sequence, G dwarfs. This section, therefore, dis- 
cusses trends of jitter with respect to AF Ca jj, 
AM V , and B—V. 

3.1. Activity 

Fig. 3 shows function of F c& jj for stars 

with AMy < 1. There is a clear increase in a' rv 
with F Ca n, such that the most active stars con- 
sistently show more than 4 ms _1 of jitter. A few 
stars with Ai 7 ' Ca jj < 1 appear to show more than 
~ 20 ms _1 of jitter; these stars may harbor planets 
whose orbital periods have not yet become appar- 
ent from the radial velocity time series. However, 
two of these stars may be very young (< 1 Gyr) 
according to Valenti & Fischer (2005), so one can- 
not discount the possibility that these devalues 
are due purely to astrophysical jitter. 

This work employs the following method to fit a 
polynomial to these data, which contain a few high 
outliers and points for which only an upper limit 
can be set. First, I divide the distribution into four 
evenly spaced bins in logAF Ca jj. Within each 
bin I then find the median value, 20th percentile, 
and 80th percentile of a' rv (represented by large 
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gray filled circles in Fig. 3), yielding three sets 
of points to which I fit second-order polynomials. 
In Figure 3, dashed lines represent the fits to the 
20th and 80th percentile points, and the solid line 
represents the avereage of these fits, representing 
the model of the median. The Appendix describes 
the definition of a median and a percentile in a 
data set containing upper limits. Note that the 
choice of a second-order polynomial here is not 
physically or statistically motivated, but is merely 
chosen to be a simple, accurate description of the 
distribution with a smooth transition to the low- 
activity subsample. 

There is very little change in the median jitter 
for stars with AF Ca jj < 0.6, so these fits are only 
applicable for stars with AF Ca jj > 0.6. Of the 72 
such stars in the sample, 53% lie below the solid 
line in Fig. 3, 80% lie below the 80th percentile fit, 
and 20% lie below the 20th percentile fit. Table 2 
contains coefficients for these polynomials. 

The result here, that jitter increases with ac- 
tivity is consistent with the result of Saar et al. 
(2003), who found a' rv oc (R / HK ) b , with b = L1 
for G and K dwarfs and b = 1.7 for F dwarfs. 
Although here the fit is parabolic and is applied 
to a different metric, one can compare the results 
by fitting a line to a' rv versus \ogR' HK for the F 
stars and the G and K stars separately, using the 
ASURV statistical analysis package rev 1.2 Laval- 
ley, Isobe, & Feigelson (1992), which implements 
the methods presented in Isobe, Feigelson, & Nel- 
son (1986) to handle the upper limits. The best-fit 
values for found b by applying the Buckley-James 
method and the parametric EM algorithm are 
b = 1.6 for F dwarfs (B-V < 0.63, AM V < 1) and 
b = 0.80 for G and K dwarfs (0.63 < B-V < 1.4, 
AMy < 1). While these results differ somewhat 
from the results of Saar et al., the strong color 
dependence remains. 

However, using the activity metric AF Ca jj in- 
stead, however, with its extra factor of T* ff , this 
color dependence is somwhat reduced: a similar 
fit to the relation a' rv oc (logAF Ca jj) c returns 
c = 0.59 for G and K dwarfs and c = 1.01 for 
F dwarfs. The choice of AF Ca jt thus reduces the 
color dependence of the jitter-activity relation, as 
shown in Fig. 4. 



3.2. Stellar Evolution 

Fig. 5 shows (T^^j . clS cl function of AMy for 
stars with AF Ca n < 0.6. Stars with AMy > 2 
have higher typical cr^-values than stars on the 
main sequence. The trend here is weak, so in this 
case, the model consists of three broad groups, 
each modeled independently: main-sequence stars 
(AM V < 1), slightly evolved stars (1 < AMy < 
1.75), and subgiants (AMy > 1.75). 

Fig. 6 shows the distribution of a' rv values for 
these three bins in AMy, with shaded regions rep- 
resenting those stars for with a' rv exceeds cr moas /2 
(that is, upper limits on a' rv ) (Eq. 7). The median, 
20th, and 80th percentile levels for these distribu- 
tions appear as vertical gray bars in Fig. 6, and 
numerically in Table 2. Again, the appendix de- 
scribes the statistics of medians for these data sets. 

Fig. 6 shows that subgiants exhibit higher a' rv 
values than their main-sequence counterparts, 
having a median a' rv value > 5ms _1 , significantly 
higher than the non-astrophysical sources of noise 
discussed in § 2.4.2. In fact, no subgiant in the 
sample performs better than 3 ms _1 . These stars 
are therefore poor velocity standards, due to jitter. 

The distribution of slightly evolved stars is very 
similar to that of main sequence stars. This work 
nonetheless treats them separately here, because 
their activity levels, normally a good indicator of 
jitter, are generally lower and therefore not com- 
parable to main-sequence stars (Wright 2004). 
Their numbers are too few for further subdivi- 
sions into bins of color or activity, so the model 
given here simply incorporates median a' rv values 
for this group as a whole. 

3.3. Color 

Fig. 7 shows the distribution of a' rv values for 
stars in three B — V ranges (corresponding roughly 
to spectral types F, G and K, and M) that have 
AMy < 1 and F Ca jj < 0.6. Here again, verti- 
cal bars show the values of the 20th percentile, 
median, and 80th percentile, and shaded regions 
represent upper limits. 

Jitter in these unevolved inactive stars is not 
a well-determined, smooth function of color, but, 
broadly speaking, three distinct B — V groups ex- 
ist. Blue stars (B — V < 0.6) generally show ~ 5 
ms _1 of jitter, while G and K dwarfs exhibit a 
median a' rv value of 3 ms _1 . The large number 
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of shaded stars in the distribution of F stars in 
Fig. 7 reflects the fact that, as noted in § 2.1, 
Planet Search stars do not receive more exposure 
time than their jitter levels warrant, to keep er- 
rors photon-limited. Note that this practice does 
not artificially inflate the estimate of jitter in these 
stars because the method of calculating medians 
accounts for these upper limits (Appendix). 

The a' rv distribution of M dwarfs exhibits a tail 
of high jitter values [a' rv ~ 7 ms _1 ), which most 
likely represents systematic errors in the deconvo- 
lution of their spectra, made difficult by the lack of 
continuum in such cool stars. The large shaded re- 
gions of the M dwarf distribution are due to these 
stars' faintness: as noted in § 2.1, no Planet Search 
star receives more than 10 minutes of exposure at 
a time, so faint M dwarfs will yield lower S /N ra- 
tios and thus have photon-limited precision. 

G and K dwarfs yield the best radial-velocity 
performance, with over 30 members exhibiting rms 
velocities below 2 ms _1 . It is fortunate that the 
best performing stars in the sample are those that 
most closely resemble the Sun. 

3.4. Other Regions 

Stars that are both active and either evolved 
or very red (B — V > 1.4) are not discussed here, 
because there are too few sample stars to properly 
determine the behavior of such stars. 

Only two sample stars meet the criteria for 
study here (§ 2.1) and have AF Ca n > 0.6 and 
B-V > 1.4. They have a' rv values of 20 and 11 
ms _1 , both higher than the values predicted by the 
median curve in Fig. 3, which are ~ 5 ms _1 . This 
indicates that red, active stars may have substan- 
tially more jitter than their bluer or more-inactive 
counterparts. 

The region of AF Ca jj-AMy space correspond- 
ing to active (AF CaII > 0.6), evolved (AM V > 1) 
stars is poorly populated, not just in this sam- 
ple, but also in other surveys of Sun like stars 
(Wright 2004). The sole such star in this sam- 
ple has <j' rv = 21 ms _1 , substantially higher than 
the 6 ms~ predicted by the u^-activity relation 
here. Such stars often represent massive early F 
stars that have entered or are entering the sub- 
giant phase before losing substantial amounts of 
angular momentum (do Nascimento et al. 2003). 



4. The Jitter Model 

Table 2 contains the a' rv values that are typ- 
ical for various regions given in [B — V)-F c jj- 
AMy space. The jitter model assumes that a 
given star will exhibit levels of jitter similar to 
those of other stars in its region. When the Planet 
Search pipeline fits Keplerian orbits to radial ve- 
locity time series, it incorporates the jitter-model 
estimate of the jitter into its estimate of the errors. 
Furthermore, the model also consists of interpret- 
ing the median values in the appropriate region 
of Table 2 as the estimated jitter, but with the 
following minor modifications. 

To prevent large discontinuities in the jitter 
model, the model linearly interpolates the me- 
dian jitter for stars within 0.1 mag of the B — V 
boundaries. For instance, an unevolved, inactive 
star with B — V = 0.65 would have predicted jit- 
ter of 3.725 ms _1 : since it is within 0.1 mag of 
the boundary between the fourth and fifth en- 
tries in Table 2, the model interpolates between 
the median value for the B — V< 0.6 and 0.6 < 
B — V < 1.4 regions. The model also similarly in- 
terpolates between the first and second regions of 
Table 2 for stars with 1.65 < AM V < 1.85. The 
model also sets the predicted jitter of active stars 
(F Ca n > 0.6) to the greater of the predicted jit- 
ter and that predicted for an inactive star with the 
same B — V and My. 

Both this jitter model and the regions and val- 
ues in Table 2 have evolved and become more 
sophisticated over time, so the model described 
here may not agree precisely with jitter estimates 
published in previous work by the California and 
Carnegie Planet Search. 

This model is generally applicable for most 
Planet Search stars. Since it was derived using 
Keck data, in particular data before the Septem- 
ber 2004 upgrade of the HIRES CCD, it accounts 
for systematic errors that may not be present in ra- 
dial velocity data from other telescopes. For jitter 
values > 3ms _1 , however, it predominantly mea- 
sures astrophysical jitter, and so is generally appli- 
cable. The presence of undetected planets around 
the sample stars means that these values may be 
slightly overestimated, but the use of medians mit- 
igates this effect. 

Jitter is more important to the detection of ex- 
oplancts than the characterization of their orbits. 
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Table 2 

Jitter values for different regions in AF CaII , B-V, AM v space. Values in curly braces 

ARE COEFFICIENTS IN A POLYNOMIAL FIT: FOR AN ENTRY {a ,ai,a 2 }, \oga' rv = J2i a i( B ~ V T ■ 





Region 






<,(m S -') 




A My 


AF CaII 


B-V 


20th percentile 


median 


80th percentile 


1.75 < A My < 4 


< 0.6 




4.6 


5.7 


9.5 


1 < A My < 1.75 


< 0.6 




2.6 


4.0 


6.2 


< 1 


6 ^ AF CaII< 3 


B-V < 1.4 


{0.64,0.92,0.55} 


{0.83,0.96,0.61} 


{1.02, 1.01,0.68} 


< 1 


< 0.6 


< 0.6 


2.1 


4.4 


6.8 


< 1 


< 0.6 


0.6 < B-V < 1.4 


2.3 


3.5 


5.3 


< 1 


< 0.6 


> 1.4 


2.5 


3.9 


6.7 



The effects of including jitter in one's noise esti- employed R' HK as a metric, 

mate on best-fit orbital parameters for a single- 
planet system are generally very small (< 1%), 
in part because the errors are typically already 
nearly equal. Jitter plays a more important role 
in determining whether radial velocity variations 
are of stellar or Keplerian origin. To the extent 
that the jitter measured here is astrophysical and 
not due to systematics or still-undetected planets, 
the 2 — 5 m s -1 of jitter exhibited by most inac- 
tive stars represents a minimum detection thresh- 
old below which planets will be very difficult to 
detect. 

5. Conclusions 

This work concerns a sample of the California 
and Carnegie Planet Search data stars represent- 
ing stars without obvious or detected planets ob- 
served at the Keck Observatory. This sample de- 
scribes the typical radial velocity stability, mea- 
sured as a' rv , of most stars in the Planet Search. 
The metric <j' rv includes the effects of undetected 
planets, systematic errors (< 3ms -1 ), and astro- 
physical jitter. 

Among inactive stars F, G and K, and M dwarfs 
show sightly different amounts of radial velocity 
jitter, from 3.5 ms -1 for G and K dwarfs to nearly 
5 ms _1 for M dwarfs. Evolved stars (AM V > 1.75) 
show significant amounts of jitter, while those only 
slightly evolved appear more consistent with main 
sequence stars. 

Active stars show markedly more jitter (5 — 10 
m s -1 ) than inactive stars, and this jitter increases 
with activity. The activity metric F Ca jj of Rut- 
ten can provide an estimate of jitter without the 
strong color dependence of previous models that 
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A. Medians and Percentiles of Data Sets That Include Limits 

The median is defined as the A/2th largest element of a sorted TV-element data set. The presence of 
measurements for which one can only set limits (nondetections, for instance) complicates matters somewhat, 
but a median can still be similarly defined. 

Consider a sorted iV-element data set x, comprising elements Xi where < i < N. Each Xi represents a 
measurement that is either exact (a data point), an upper limit, or a lower limit. The set x should be sorted 
such that for otherwise equal elements, an upper limit is less than a data point, which in turn is less than a 
lower limit. (Note that these definitions apply only to the sorting procedure; different comparison conditions 
apply below.) One seeks to test the value of each element Xi as a candidate median by comparing a,, the 
number of elements unambiguously above the value of x i} to bi, the number unambiguously below it. 

For the purposes of computing 6,, element Xj is counted as being below the value of Xi if, in the sorted 
array: (1) Xj is a data point and j < i or (2) xj is an upper limit and j < i. Note that when Xi represents an 
upper limit, Xi is less than the value of Xi by definition, and so contributes to bi, and that when xj represents 
a lower limit, it does not contribute to bi at all. Similarly, Xj contributes to a, if j > i when Xj is a data 
point, or if j > i when Xj IS cl lower limit. 

The median is then is the value of the ith element of x such that a* is most nearly equal to bi. If more 
than one such element exists, these can simply be averaged. In some cases, such as a data set consisting 
entirely of upper limits, no median can be defined. Similarly, the pth percentile can be defined as the value 
of Xi such that bi/(ai + bi) is most nearly p/100. Naturally, the 50th percentile is equivalent to the median. 

For example, consider the appropriately sorted data set x = {< 1,5,5, < 7,10,50}, in which the first 
and fourth elements represent upper limits. The values of a for this set are {4,3,2,2,1,0}, and similarly 
b = {1, 1, 2, 4, 4, 5}. Thus, the median of this set is 5 (= X2, since 02 = 62), the 20th percentile is 1 (= xo, since 
bo/{ao+bo) = 1/5 = p/100 whenp = 20), and the 80th percentile is 10 (= X4, since 64/ (04+64) = 4/5 = p/100 
when p = 80). 

I am indebted to Geoff Marcy and the California and Carnegie Planet Search for the data, guidance, 
and advice that made this work possible. I also thank the anonymous referee for the extremely timely and 
constructive comments that improved this paper. 

I wish to recognize and acknowledge the very significant cultural role and reverence that the summit of 
Mauna Kea has always had within the indigenous Hawaiian community. We are most fortunate to have the 
opportunity to conduct observations from this mountain. 
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NASA's Astrophysics Data System Bibliographic Services. It was made possible by the generous support of 
Sun Microsystems, the NSF, and NASA through grant NAG5-12090. 
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Fig. 3. — Activity, AF Ca n, vs. the observed jitter, cr' rv for stars with AMy < 1. Typical errors on data 
points are around 1.5 ms^ 1 . Upper limits are represented by arrows. Large gray disks represent the 20th 
percentile, median, and 80th percentile in each of 4 equally spaced bins in log A_F Ca jj from 0.3 < AF Ca jj < 3. 
The dashed lines represent the parabolic fits to these points in the range 0.6 < AF Ca tt < 3, where the model 
is applicable. 
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Fig. 4. — Jitter as a function of activity using two different activity metrics. Solid lines represent fits using 
only F stars, dashed lines those using only the G and K stars. F stars show a steeper activity dependence 
than G and K stars with logi?[j K as a metric (left panel), but the difference is smaller when F Ca j] is the 
metric (right panel). These fits are provided for comparison only, and arc not used in the jitter model 
presented here. 
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Fig. 5. — Evolution, measured as AMy vs. cr^, for stars with F Ca n < 0.6. Typical errors on data points 
are around 1.5 ms _1 . Arrows represent upper limits. The most evolved stars exhibit large amounts of jitter. 
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Fig. 6. — Distribution of a' rv , the observed jitter, for inactive (A_F Ca jj < 0.6) stars in bins of AMy. The 
top panel includes only stars for which 0.6 < B — V < 1.4, matching the majority of stars in the bottom two 
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Fig. 7. — Distribution of a' rv , the observed jitter, for inactive (AF Ca jj < 0.6), unevolved (AMy < 1) stars 
in three color bins. The intermediate color stars show the lowest median jitter. Shaded regions indicate stars 
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